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The magnetic energy of nonlinear Alfven waves in com- 
pressible plasmas may be ponderomotively coupled only to 
ion-acoustic quasi-modes which modulate the wave phase ve- 
locity and cause wave-front steepening. In the collisionless 
plasma with /3 ^ 0, the dynamics of nonlinear Alfven wave is 
also affected by the resonant particle- wave interactions. Upon 
relatively rapid evolution (compared to the particle bounce 
time) , the quasi-stationary wave structures, identical to the so 
called (Alfvenic) Rotational Discontinuities, form, the emer- 
gence and dynamics of which has not been previously under- 
stood. Collisionless (Landau) dissipation of nonlinear Alfven 
waves is also a plausible and natural mechanism of the solar 
wind heating. Considering a strong, compressible, Alfvenic 
turbulence as an ensemble of randomly interacting Alfvenic 
discontinuities and nonlinear waves, it is shown that there ex- 
ist two distinct phases of turbulence. What phase realizes de- 
pends on whether this collisionless damping is strong enough 
to provide adequate energy sink at all scales and, thus, to sup- 
port a steady-state cascade of the wave energy. In long-time 
asymptotics, however, the particle distribution function is af- 
fected by the wave magnetic fields. In this regime of nonlinear 
Landau damping, resonant particles are trapped in the quasi- 
stationary Alfvenic discontinuities, giving rise to a formation 
of a plateau on the distribution function and quenching col- 
lisionless damping. Using the virial theorem for trapped par- 
ticles, it is analytically demonstrated that their effect on the 
nonlinear dynamics of such discontinuities is non-trivial and 
forces a significant departure of the theory from the conven- 
tional paradigm. 

PACS numbers: 52.35. Mw, 52.35.Ra, 96.50.Bh, 96.50.Ci 



I. INTRODUCTION 

Alfven waves are the most 'robust' plasma oscillations 
in magnetized systems. They are encountered in a wide 
variety of astrophysical objects such as stellar and solar 
winds, interstellar medium, hot accretion flows, magne- 
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tospheres of planets and magnetized stars, etc.. It is in- 
teresting that in incompressible plasma, V • v = 0, fluid 
and magnetic nonlinearities in a finite-amplitude Alfven 
wave cancel each other exactly, so that it obeys the same 
dispersion relation as the linear wave of infinitesimal am- 
plitude. 

Such a mutual cancellation of Reynolds and magnetic 
stresses in the wave breaks down when plasma is com- 
pressible. Indeed, the magnetic field pressure, B^/Stt, 
associated with the wave field ponderomotively modu- 
lates the local plasma density, n, which, in turn, affects 
the phase velocity of the wave. This process may also be 
viewed as an effective parametric coupling of an Alfven 
wave to acoustic waves in plasma. Such acoustic modes, 
however, are not plasma eigenmodes. They are driven by 
the Alfven wave and, hence, propagate with the Alfven 
speed, which is, in general, different from the sound speed 
in plasma. The nonlinear coupling of the phase speed of 
the wave to its amplitude in the nonlinear Alfven waves 
is ultimately responsible for the wave-front steepening 
and formation of collisionless shocks in astrophysical and 
space plasmas. 

In studies of nonlinear waves it is customary to assume 
'weak nonlinearity', i.e., the wave amplitude is finite but 
still small compared to the ambient field. This means 
that a wave is linear in the leading order of expansion in 
wave amplitude. All nonlinear effects appear as a grad- 
ual, large-scale modulation of the wave amplitude. In 
other words, one studies slow evolution of an envelope of 
a wave packet. 

It is well known that in collisionless, incompressible 
plasma, low-frequency Alfven waves are undamped or 
very weakly damped since collisionless dissipation is in- 
efficient. Ion-cyclotron damping is always weak for low- 
frequency waves and Landau damping, which could be 
important, is absent because there are no longitudinal 
electric and magnetic field perturbations in the Alfven 
wave. Equivalently, an Alfven wave is strictly transverse 
in the drift approximation and thus cannot affect the par- 
ticle motions along the field lines. (Note that weak Lan- 
dau dissipation of obliquely propagating Alfven waves 
occurs on scales comparable to the ion Larmor radius. 
Pi, where the drift approximation breaks down. This 
damping is really weak since typically k±pi <C 1.) 

In contrast, a nonlinear Alfven wave in compressible 
plasma is nonlinearly coupled to acoustic-type oscilla- 
tions which are accompanied (in a two-fluid model) by 
longitudinal electric fields and are damped via the Lan- 
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dau mechanism. Hence, dissipation enters the dynam- 
ics of nonUnear Alfven waves nonUnearly, even in the 
regime when the resonant particle response is calculated 
in the linear approximation, i.e., neglecting trapping of 
resonant particles. Another important and quite unusual 
feature is the nonlocality of the damping. It is associ- 
ated with the finite transit time of a resonant particle 
through the wave envelope modulation and represented 
by an integral operator in the wave evolution equation. 

Unusual properties of dissipation in the nonlinear 
Alfven waves naturally result in the following peculiar 
features of these waves. First of all, it was shown that 
Alfven waves do not always damp to zero. Instead, 
they evolve into the localized, coherent, quasi- stationary 
structures referred to as directional and rotational discon- 
tinuities. These are the narrow regions where the wave 
magnetic field vector rapidly rotates through some an- 
gle and, respectively, are or are not accompanied by the 
wave amplitude variations. These discontinuities con- 
nect the parts of a wave packet of different polarizations 
or wave phases. Second, among the Alfvenic discontinu- 
ities, there are solutions which constitute a new class of 
shock waves, — the collisionless dissipative shocks. Such 
shocks form when nonlinear wave steepening is balanced 
by collisionless damping process (and not by the wave 
dispersion, as in usual collisionless shocks). Due to scale 
invariance of Landau damping, the shock width is de- 
termined by dispersion, however, which is the only char- 
acteristic scale in the problem. It should be emphasized 
that Alfvenic rotational discontinuities have been discov- 
ered by in situ measurements |^-|^ of magnetic fields in 
the solar wind. Despite intense observational and theo- 
retical efforts these years, a comprehensive quantitative 
theory which explains the emergence and dynamics of 
Alfvenic discontinuities was proposed only recently [^-^ . 
All the results are in excellent agreement with the ob- 
servational data and numerical hybrid simulations [0,1). 
Third, preliminary study of a one-dimensional model of 
strong, compressible Alfvenic turbulence suggests the 
existence of two distinct phases or states of the turbu- 
lence, depending on how strong the collisionless damping 
is. At last, resonant particles, which were responsible for 
the wave damping in the linear Landau process, appear 
to be trapped in the wave potential at large times and 
the linear theory of Landau damping is invalid. Qual- 
itatively, the damping rate oscillates with roughly the 
bounce frequency and gradually decreases to zero un- 
til resonant particles phase mix and form a plateau on 
the distribution function. A recently proposed asymp- 
totic theory of particle trapping |10| predicts non-trivial 
nonlinear contributions of resonant particles to nonlinear 
Alfven wave dynamics. They are the following. (1) The 
power of the KNLS nonlinearity associated with resonant 
particles increases to fourth order when trapped parti- 
cles are accounted for. (2) The effective coupling is now 
proportional to the curvature of the particle distribution 



function (PDF) at the resonant velocity, /g (wa), and not 
to its slope, /^{va), as in linear theory. (3) The phase- 
space density of trapped particles is controlled by plasma 
/3. Thus, the structure of quasi-stationary Alfvenic soli- 
tons and discontinuities may be modified in r — *■ (X) limit. 

The rest of this paper is organized as follows. In Sec. 
H a semi-qualitative derivation of the envelope evolution 
equation for finite- amplitude Alfven waves, referred to 
as the Derivative Nonlinear Schrodinger (DNLS) equa- 
tion, wave is presented. In Sec. Ill we generalize it to 



include the resonant particles effect (Landau damping) 
and obtain the kinetically modified derivative nonlinear 
Schrodinger equation (KNLS). We discuss numerical so- 
lutions for this equation which appear to be the Alfvenic 
Rotational Discontinuities that observed in situ in the 
solar wind. We also discuss the energy dissipation time 
scales for various types of nonlinear Alfven waves. In Sec. 
IV the noisy-KNLS model of strong, compressible magne- 
tohydrodynamic (MHD) Alfvenic turbulence is discussed. 
In Sec. a self-consistent treatment of the nonpertur- 
bative effects of particle response (particle trapping) on 
the nonlinear Alfven wave dynamics is presented. 



II. NONLINEAR ALFVEN WAVES: THE DNLS 
EQUATION 

The usual derivation of the evolution equation of the 
nonlinear Alfven waves is based on the multiple time scale 
expansion. It assumes that the envelope of an Alfven 
wave varies on much longer time-scale, r = {Bq / B±)^t ^ 
t, than the linear wave period, t = lu^^ = (k^^VA)^^, 
where Bq and B are the unperturbed magnetic field and 
its perturbations, va = Bq / ^ Airmin is the Alfven speed, 
rui is the ion mass, and || and _L designate components 
parallel and perpendicular to the ambient (unperturbed) 
magnetic field. Such a derivation has been performed by 
different authors and we refer the reader to the original 
literature ||l^,|l^. To provide physical insights, we give 
here a transparent, semi-qualitative derivation instead. 

A wave equation of a small amplitude Alfven wave may 
be written as 



{LJ~k\iVA [l±k\\VA/2ni\) Bj 



0, 



(1) 



where fli — eBa/niiC is the ion-cyclotron frequency. The 
last term is simply weak dispersion, fc||t;A/r2,j ^ 1, due 
to a finite Larmor radius. A ponderomotive force creates 
density fluctuations, 5n, so that the Alfven speed is no 
longer constant: 



Bq / 1 5n 

VA — — 1^^^^=^^^^= — VAQ 1 

A'Kmi{nQ + 5n) \ '2 no 



(2) 



We assume that the wave propagates along the magnetic 
field in z-direction (to the right) and its amplitude is 
changing gradually, B± = B±{z — VAot^r). Going back 
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to real-space representation, this suggests the following 
replacement: 
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We can neglect the variation of va in the small dispersion 
term as a next order effect. From Eq. (|l|) it is easy to 
write: 



2 dz \no 



0. 



(4) 



Hereafter we omit subscript "0" in va- The perturbed 
magnetic field here is a complex quantity and is defined 
as B± = Bx±iBy. The density perturbations are caused 
by the ponderomotive force exerted on plasma by the 
fluctuating magnetic field of the wave. They arc obtained 
from the continuity, momentum balance equations, and 
equation of state which yield the equation for ion-acoustic 
waves with source: 



a<2 



' dz' 



Sn 
no 



2 dz 
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where — ^{T^. + Ti) jmi is the ion- acoustic speed, 7 = 3 
is the polytropic constant, and p = n{Te + Ti) and rrie <C 
nii were used in derivation. Since, ion-acoustic waves are 
driven by the Alfven wave, the density perturbation is a 
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and Eq. (H) yields: 
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where P = c^/v^ is roughly the ratio of plasma pressure 
to magnetic pressure. Substituting the density perturba- 
tion into the equation for the wave magnetic field (jj) and 
defining b = B±/Bo, we obtain the DNLS equation: 
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This equation describes evolution of planar nonlinear 
Alfven waves propagating in one dimension along the 
ambient magnetic field. It was shown that this equa- 
tion is solvable [l3| and admits soliton and cnoidal wave 
solutions. 

A numerical solution of the initial value problem for 
the DNLS is shown in Fig. |l|a. It is natural to assume 
that nonlinear waves emerge from small-amplitude (lin- 
ear) ones. Thus the most general class of initial pro- 
files are the finite-amplitude, periodic waves of different 
helicities (polarizations). In the case of Fig. [^a, an ini- 
tially linearly polarized, sinusoidal initial wave profile has 
been chosen. (The range of dimensionless coordinate, 
< < 1024 corresponds to the number harmonics used 
in calculation, —512 < k < 512.) One can see that wave 
front indeed steepens at early times, r ~ 2 ~ 400ri~^. 



Later, at t '--^ 5, when the width of the front becomes 
of order fli/vA dispersion limits nonlinear steepening 
and produces small-scale, oscillatory wave- like structures 
with significant small-scale energy component. At later 
times, the magnetic field is highly irregular indicating 
strong nonlinear Alfven wave turbulence. Note that no 
(quasi- stationary) Alfvenic discontinuities emerge. 



III. NONLINEAR ALFVEN WAVES WITH 
LANDAU-TYPE DISSIPATION: THE KNLS 
EQUATION 

It is known that kinetic effects like Landau damp- 
ing are absent from the MHD approximation. To in- 
clude them into the nonlinear Alfven wave model self- 
consistently, a fully kinetic treatment is needed jlj] . The 
results are often too complicated and obscuring, however. 
There are two ways to go around the difficulties. First 
way is to use a hybrid MHD-kinetic approach in which 
one performs the ad hoc calculations the quantities which 
are known to be affected by kinetics (i.e., the density 
perturbations, Sn) and plugs them into the MHD model 
|l5|"p^. Second way is to make some self-consistent clo- 
sures in deriving the MHD equations to modify them so 
that to mimic collisionless processes |18 . Then these new 
MHD equations may be used to re-derive the nonlinear 
Alfven wave equation Here we again give a simple, 
semi-qualitative derivation. For rigorous analysis, we re- 
fer reader to the original literature ||5|,p"5|-p^ . 

To include kinetic effects, we need to modify the equa- 
tion for ion-acoustic waves only, while Eq. (^) remains 
unchanged. The equation of damped ion-acoustic modes 
may be formally written as 
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d 
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where 7 is a damping rate. If the damping rate is a func- 
tion of k in Fourier space, 7 in real space is represented 
by an operator acting on a wave field. We now recall 
that in /c-space Landau damping rate is proportional to 
l/|fc|[|. Thus in real space we write 



7 = VAX>i 



d_ 

dz 



(9) 



where 7i[/] = ^ p3^/(z')dz' is the integral Hilbert 
operator acting on a function / [V denotes principal 
value integration), which is equal to (ifc||/|A:|| |)/fc in 
Fourier space and, thus, ensures correct /c-scaling of col- 
lisionless damping. Here x is a parameter which contains 
all insignificant physics such as the structure of a PDF, 
etc., and depends on the closure chosen. In the model 
of Ref . Q , X has the physical meaning of a thermocon- 
ductivity coefficient. Again, we assume a traveling wave 
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solution for the density perturbations, 5n = Sn{z — VAt). 
Than we have 



no 



5n 
no 



(10) 



To solve this equation we use the identity TCH = —1, 
which has a physical meaning of time-reversibility of Lan- 
dau damping [||. Acting on Eq. ( p^ with H, excluding 
7i[5n/no] and substituting into Eq. we obtain the 
KNLS equation: 



db 5 / , 
^ + '"A^ [mil 



m2bH[\b\'' 



-n 

(11) 



where mi = (l-/3)/4A, ma = -x/4A, A = (l-/3)^+x^ 
are coefficients, which are mo del- dependent. The rigor- 
ous derivation j|] yields 



mi 



m2 



1 (l-/3*)+x'(l-p*/7) 
/3*)2+x'(l-/3*/7)'^ 
X(3*{j-l)/j 



4(1 
1 

"4(l-/3*)2 + x2(l-/3V7)^ 



(12a) 
(12b) 



where = (Te + r,)/3/T, ~ (Te/T,)/3 for T, < Te, and x 
is the best fit parameter of the closure to the exact kinetic 
particle response [|l8| with a temperature correction, x = 
y8^Mexp{-l//3}(re/T,)3/2 exp{(T, - Te)/2T,}. The 
dependence of mi and m2 vs. 1//3 is drawn in Fig. |^. 
Negative sign of m2 indicates damping. The damping of 
an Alfven wave is strongest in warm, isothermal plasmas, 
/3 - 1, Te ~ T,. The KNLS equation straightforwardly 
generalizes to the case of obliquely propagating waves. 
In fact, Eq. ([Ill) remains unchanged, but the wave field 



is re-defined as follows: b = {Bj. + iBy + BosinQ)/Bo, 
where Q is the obliquity angle, s'mO oc k • Bq. 

Clearly, mi and m2 are functions of Te/Ti and /3. It 
can be shown that if either T^/Ti ^ 1, or /3 ^ 1, or 
/3 <C 1, the Landau damping of nonlinear Alfven waves 
becomes very weak, m2 0. Then the DNLS nonlin- 
earity dominates and mi^(l — /3)~^,i.e., changes sign 
at /3 ~ 1. In this regime, other dissiaption mechanisms 
(e.g., ion-cyclotron damping of the linear Alfven wave 
[carrier]) may become important, see Sec. [IIB. 



A. Alfvenic discontinuities 



The KNLS equation (|ll|) reduces to the DNLS for 
/3 — > 0. CoUisionless damping of nonlinear Alfven waves 
vanishes in this case. The wave dynamics is thus decou- 
pled from resonant particle effects and Alfvenic discon- 
tinuities do not emerge (unless one artificially plugs in 
diffusion, finite plasma conduction, or other coUisional 
effects, as done in some simulations). For finite /3's, the 



nonlinear Alfven wave evolution is affected by resonant 
particles and is described by the second nonlinear, inte- 
gral term in Eq. (pi]). Unlike the DNLS equation, the 
KNLS equation is not integrable because of dissipation. 
Numerical solution of this integro-differential equation 
for the same initial conditions as in Fig. 0a and = 1, 
reveals different wave dynamics, as shown in Fig. 
Instead of irregular, fluctuating wave fields, localized, 
quasi-stationary waveforms (Alfvenic discontinuities) are 
seen to form very rapidly, within r ^ 5 = lO'^ri"^. There 
are three parameters (in addition to P and Te/Ti) which 
control the wave dynamics as well as specify the type of 
discontinuity which results. These are (i) the wave helic- 
ity (i.e., the polarization type), (ii) the obliquity angle 0, 
and the angle between the polarization plane (for nearly 
linear polarizations) and k-Bg plane. 

There are three general types of quasi- stationary 
Alfvenic discontinuities, discriminated by the obliquity 
and phase jump (i.e., angle though which the magnetic 
field vector rotates). First, there is a wide class of the arc- 
type rotational discontinuities which emerge in oblique 
propagation, Q > 10°, independent of the wave helic- 
ity (polarization) . They are characterized by an arc-type 
shape diagram in the hodograph (the bx-by diagram), 
as shown in the "snap-shot" , Fig. ||. The wave phase 
jump through the discontinuity is Afjj < tt. Second, there 
is a class of the S-type directional discontinuities which 
emerge in quasi-parallel propagation, 0° , only from 
linearly and almost linearly polarized initial perturba- 
tions. They are characterized by a remarkable S-shaped 
hodograph, as shown in Fig. 0. They are called direc- 
tional (not rotational) discontinuities because the phase 
jump is accompanied by moderate amplitude variation. 
The wave phase changes by A(j> = ir through the the dis- 
continuity. Third, there is a narrow class of arc-type ro- 
tational discontinuities propagating parallel to the mag- 
netic field. They emerge from waves with elliptical polar- 
ization only. The wave phase change is A(j) = tt. Parallel 
propagating circularly polarized waves (helicity equals to 
unity) do not evolve to discontinuities and are decoupled 
from dissipation. 



B. Comparison with observations 

There are two regimes of the Alfven wave evolution, 
depending on the values of Te/Ti and (3. Although there 
is no sharp boundary between the regimes, a qualita- 
tive insight may be gained from Fig. ^, where the curve 
|mi| ~ |m2| is plotted. Kinetic damping is important for 
Te/Ti and /3 from the region labeled "hydrodynamic" . 
There |mi| < |m2| and Alfvenic rotational discontinu- 
ities rapidly form. We emphasize the remarkably good 
correspondence of the KNLS solutions and the solar wind 
observational data ||l|-^. Recent 2^ hybrid code simula- 
tions are also in excellent agreement with the results of 
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the KNLS theory and thus support the idea that 
the mechanism of formation of Alfvenic discontinuities is 
the combined effect of resonant wave-particle interactions 
and nonhnear wave steepening. 

In the opposite case, |mi| >> |m2| (labeled "bursty" 
in Fig. 1^), the Landau damping is weak and may be 
neglected compared to the ion-cyclotron one. Nonlin- 
ear wave steepening cascades wave energy to the scales 
k ^ max{\mi\,\m2\}\b\'^fli/vA, as estimated from Eq. 
(pi]). Thus, the number of particles being in the cy- 
clotron resonance, k{vA — v) — fli = 0, may be large 
if 7712 — > 0, /3 ~ 1. Nonlinear dynamics in this case is 
governed by the so called DNLS-Burgers equation [p^ . 
This theory explains well the emergence and evolution 
of Alfven shock trains (also referred to as "shocklets") 
|pO| , which were observed upstream of the bow-shocks of 
planets and comets (see e.g., Ref. [^). 

C. Energy dissipation 

Resonant particle effects make the otherwise almost 
non-dissipative Alfven waves to damp. Fig. |^ shows tem- 
poral evolution of magnetic energy associated with dif- 
ferent structures (Alfvenic discontinuities). Clearly, the 
energy dissipation rate in the case of quasi-parallel propa- 
gation is sensitive to the initial polarization, i.e., helicity, 
which is, roughly, a measure of asymmetry of a spec- 
trum between +k and —k harmonics. Crude estimates 
1^ of the typical damping time are: (i) for linearly polar- 
ized waves (zero helicity), e.g., the S-type discontinuities, 
Tiin cx |to2|~^ and the damping is algebraic, |6p cx t~^/^, 
rather than exponential, (ii) for elliptically polarized 
waves Toll cx |to2As|"\ where ~ X] (l^+feP - l^-feP) 
is a measure of spectrum asymmetry, and (iii) for circu- 
larly polarized waves (helicity unity) Tcir — > cx) as they 
do not damp. In the case of oblique propagation, the 
ambient field enters the damping rate and for all wave 
helicities Tobi cx |r7i2 sin^ 9|^^. From Fig. || one sees 
that the a typical time-scale of the wave damping is 
T ^ 100 ~ 10'*O^^, which is (roughly) a few hours for 
the solar wind conditions. 

Given the spectrum of magnetic field energy pertur- 
bations, Ek = (|&P)fe, the instantaneous rate of energy 
change is calculated exactly |^J2^ from Eq. ( pT|) to yield: 

— =m2/ |i?fc||nfc|||dfc|| . (13) 

Note that this equation is equally applicable to the sin- 
gle wave case as well as to the compressible Alfvenic tur- 
bulence in general, provided the turbulence spectrum is 
known. 

Typical dissipation timescale is several hours for the 
solar wind conditions. Hence, coUisionless dissipation of 



nonlinear Alfven waves a plausible (and natural) mech- 
anism of the solar wind heating within one astronomical 
unit. 

IV. COMPRESSIBLE TURBULENCE OF 
NONLINEAR ALFVEN WAVES 

There are two approaches in the studies of turbulence, 
cascade and structure-based. The cascade-type theories 
usually treat turbulence as a "soup of eddies" , — a col- 
lection of plasma excitations with random phases which 
interact with each other (usually on matching scales) and 
produce a cascade of energy. The structure-based-type 
approach, one of the examples of which is the noisy Burg- 
ers model in hydrodynamics, is based on the study of 
nonlinear evolution equations with external noise drive. 
It thus studies the turbulence of coherent structures, such 
as shocks, solitons, nonlinear (cnoidal) waves, etc., gen- 
erated and interacting with each other randomly. Such a 
process may be referred to as "coherent cascade" , mean- 
ing that modes (eddies) at all scales interact coherently 
(in phase) and, thus, experience stronger interactions due 
to the cumulative effect. Its is known that the first ap- 
proach is more suitable for studies of weak turbulence 
and the second one better describes the strong, highly 
nonlinear turbulence. 

The analytical study of the structure-based noisy- 
KNLS model as a generic model of coUisionless, large- 
amplitude, compressible MHD (Alfvenic) turbulence is 
presented in Ref. Stationarity is maintained via 

the balance of noise and dissipative nonlinearity. The 
Fourier-transformed KNLS equation (|Tl| ) with noise 
source reads 

(—iuj + ivok + i^J,ok^) bk + iAfc^ b^' bk" b k-k'-k" 

k' .k" ^ ^ ^ 

X [mi -f im2Sgn(fc — k')] — fk , (14) 

where vq and fio are the bare phase velocity (it is nec- 
essary for a self-consistent renormalization analysis) and 
dispersion, A = 1 is the standard perturbation param- 
eter, and the function sgn(a;) = a;/|x|. Here we omit- 
ted the subscript || by k. The one loop renormaliza- 
tion group analysis of the problem with zero mean, i5- 
correlated noise with white fc-spectrum yields complex- 
valued renormalized coefficients i^turb, Mturb- The 'Re' 
and Tm' parts of ^turb are simply turbulent dispersion 
and turbulent viscosity. The real part of Wturb may be 
interpreted as wave momentum loss via interactions with 
resonant particles and its imaginary part manifests fast, 
exponential damping due to phase mixing. Further anal- 
ysis shows the existence of two different phases of turbu- 
lence. The bifurcation occurs at the point 

\m,i/m,2\ ~ 1. (15) 
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If \mi/m2\ < 1, Landau dissipation is efficient to sink 
all the injected energy during the "coherent cascade" so 
that the hydrodynamic (tj ^ 0, fc — s- 0) regime with 
no sharp fronts realizes. In this regime, the strongly in- 
teracting Alfvenic discontinuities dominate in the tur- 
bulence. In the opposite case, nonlinearity overcomes 
damping and no stationary, hydrodynamic turbulence is 
predicted. Small-scale, bursty turbulence is expected in 
this regime. Given the weak cyclotron damping, such a 
regime corresponds to the Alfvenic shocklet turbulence. 
Recalling that mi and m2 are functions of /3 and Te/T^, 
a "diagram of state" can be drawn, as shown in Fig. O. 



V. ASYMPTOTIC, r ^ oo, DYNAMICS OF 
DISSIPATIVE NONLINEAR ALFVEN WAVES 



The linear Landau damping theory used in derivation 
of the KNLS assumes a time-independent (Maxwellian) 
PDF. It is clear that, for a finite-amplitude wave, parti- 
cles which are near resonance with the wave, v va, will 
be trapped by the wave potential U{z) = B\/^imQ be- 
cause their kinetic energy (measured in the wave frame) , 
im(z) — wa)^j is less than the potential barrier, |C/,„| — 
max\U{z)\^ as in Fig. |^. Such particles experience re- 
flections at turning points zi and Z2 exchanging en- 
ergy with the wave and significantly modify the PDF 
near resonance. The damping rate of a wave oscil- 
lates in time with gradually decreasing amplitude un- 
til phase mixing results in flattening of the PDF (for 
resonant velocities) and formation of a plateau] then 
the damping rate vanishes p2[ . Thus, the linear cal- 
culation of Landau dissipation, while correct for times 
short compared to the typical bounce (trapping) time, 
T ^ Ttr — (fc|| yJU/mi)~^ ~ (fc||Uyi|6|)~-^, fails for quasi- 
stationary waveforms (Alfvenic discontinuities) on times 
T > tnl > Ttr [t^I ~ mik\\VA{B'^/BQ) is the typical 
nonlinear wave proflle evolution time]. Hence, Landau 
dissipation should be calculated non-perturbatively to de- 
termine the resonant particle response to the nonlinear 
wave. 

The nonlinear Landau damping problem is, in general, 
not analytically tractable, as it requires explicit expres- 
sions for all particle trajectories as a function of initial 
position and time. Such trajectories cannot be explicitly 
calculated for a potential of arbitrary shape. Usually, a 
full particle simulation is required. In some cases, it is 
useful to approximate the wave profile shape by a sim- 
ple analytic expression which may be assumed to per- 
sist, while the wave amplitude varies, to calculate the 
trajectories |2^. This is not the case for our problem 
because nonlinear Landau damping controls the profile 
of the Alfven wave. 

A nonperturbative, self-consistent theory of nonlinear 
wave-particle interactions was constructed in the asymp- 
totic limit, T ^ c», applying the virial theorem to deter- 



mine the PDF. The generalized KNLS again is Eq. (Q) 
with Sn = SnNu+Snu, where Stinu and Stir cx x)^[U{z)] 
are the non-resonant (bulk) and resonant particle re- 
sponses, respectively. Here /C is a new kinetic operator 
(which replaces H) acting of the wave field which con- 
tains all the information about resonant (trapped) par- 
ticle trajectories. It is interesting that the very possi- 
bility of writing the generalized KNLS equation in this 
form relies on the intrinsic time reversibility of the Vlasov 
equation, linear or nonlinear: ICK. ~ —1, see Ref. 

The resonant particle response is calculated using Li- 
ouville 's theorem ( "the PDF is constant along particle 
trajectories"): 



f{v,z,t)^fo{v± (£,4)), 



(16) 



where z^ — z^ (z, t, E; U{z)) is the initial coordinate of 
a particle of total energy E which at time t is at the 
point z and has a velocity v± {E, z). By definition Sur = 

/A.^di; so that 



6n} 



E 

(±) 



fo{v±{E,z^))-Mv±{E,z)) 
AB™, y/2m, [E - U{z)) 



dE. (17) 



Here the sum is over particles moving to the right (+) and 
to the left (-), as in Fig. |^. The integration is over the 
resonant (negative) energies of trapped particles, Um < 
E < Q with Um being the amplitude of the potential. In 
the short-time limit, r — )■ oo), Eq. (17) reduces to the 
KNLS case with IC = H and x — ''^v\fQ{vA) /rniUQ. 

To treat the t —t oo limit we employ the virial theorem^ 
which relates the average kinetic and potential energies of 
particles trapped in an adiabatically changing potential: 



2{K{z))^n{U{z)). 



(18) 



Here U{z) = U{z) — Um is a homogeneous function of 
its argument of order n, i.e., U{az) — a"U{z). Direct 
integration of Eq. (p^ yields 



{Stir) 



foivA)J^^/WW\ 



n + 2 



X {\Um\-\U{z)\)- 



3{n + 2 



■\U{z)\ 



(19) 



Note that the term cx /o(w^) vanishes identically, hence 
the damping is absent. Since {JC){JC) ^ ICK — —1, 
we can only estimate the coupling constant to be x ~ 
foi.'^A) \y\/no\ \/2/mf . The index n is formally not de- 
fined for an arbitrary potential, but it can be estimated 
numerically, given the wave profile |l^ . 

Clearly, in the asymptotic limit, t 3> Ttr, the damp- 
ing rate vanishes due to phase mixing. Nevertheless, the 
resonant particles still contribute the wave dynamics, in 
that 
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{Sur) ~ fHivAW, (20) 

thus determining a new nonlinear wave equation. 

Simple Bernstein-Green-Kruskal type analysis shows 
that the 'height' of the asymptotic plateau on the PDF 
depends on the plasma parameters through mi and 7712 
. The difference between the plateau and the initial, 
Maxwellian PDF evaluated at va is 

/plat - /Maxw|„^ ~ ^ ' 

Thus we conclude that there is under-population of 
trapped particles (void on the PDF, /piat < /Maxw) at 
low /3's and over-population (bump, /piat > /Maxw) at 
high /3's. 

Finally, assuming weak damping of a nonlinear Alfven 
wave (e.g., ion-cyclotron damping), one can see ||ic| ] 
from the conservation of the parallel adiabatic invariant 
J — § P||dz that trapped particles will condense near the 
bottom of the wave potential: Aw|| oc 1 6 1, forming a clump 
on the PDF at ti ~ 7;^- 
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FIG. 4. Same as Fig. |3| for the S-type directional disconti- 
nuity propagating parallel to the field. 



FIG. 5. Wave energy evolution. 

FIG. 6. The /J-Te/Ti-diagram of state. The region inside 
the curve corresponds to highly damped turbulence. No steep 
fronts appear. There are steep fronts in the region outside the 
curve. 



FIG. 7. Trapping potential. 
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NOTE TO EDITOR: 
ALL FIGURES SHOULD BE SCALED TO THE SIN- 
GLE COLUMN FORMAT. 
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